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Entanglement Entropy on the Cayley Tree 
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The properties of the entanglement entropy (EE) of a clean Cayley tree (CT) are studied. The 
EE shows a completely different behaviour depending on the way the CT is partitioned into two 
regions and whether we consider the ground-state or highly excited many-particle wave function. 

The ground-state EE increases logarithmically as function of number of generation if a single branch 
is pruned off the tree, while it grows exponentially if the region around the root is trimmed. On the 
other hand, in both cases the highly excited states’ EE grows exponentially. Implications of these 
results to general graphs and disordered systems are shortly discussed. 
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I. INTRODUCTION 

In recent years there has been a renewed intereslf^ in 
the problem of Anderson localization on the Cayley tree 
(CT) (in the infinite limit known as a Bethe lattice). This 
interest is mainly motivated by the connection between 
the many body localization phenomenoiP and the CT. 
Many body localization may be viewed as a localization 
problem in Fock space, where the coupling between states 
due to electron-electron interactions resembles the ctI^. 

In the context of the Anderson localization on the CT^^, 
one can envisage that the entanglement entropy (EE) will 
play an important role in clarifying the CT localization 
properties. In this EE could join other methods such as 
level spacing statistics for the CT system^^^^H^, and for 
many-body interacting systems believed to map on an 
effective CT in Fock spac^^®^. Prior to addressing the 
challenging question of the EE on a disordered one must 
clarify the behavior of the EE for a clean CT. This is our 
main goal in this paper. 

The EE is a measure of the entanglement between two 
regions, A and B, of a system which is in some pure state. 
This measure is given by the von-Neumann entropy of the 
reduced density matrix, pA^ of region A: 

5A = -Tr(p^lnp^). (1) 

It has been shown that the EE is a very useful measure to 
locate and analyze quantum phase transit ion^^^H^. The 
ground state EE typically scales like the boundary area 
of the regiorP^^. Thus, for one dimensional systems, 
one would expect that the EE will be constant. This is 
correct for gaped systems, but for metallic systems there 
is a logarithmic correction, and EE grows like ln(I/^), 
where La is the size of the sub-system A. For insulators 
the EE should not depend on the region’s size. Therefore 
for disordered systems the EE can give an indication for 
the localization length we expect the EE to grow log¬ 
arithmically for ^ ^ and to saturate for ^ ^ I/a • 

Applying this picture to the CT is not straight for¬ 
ward. Unlike the situation for regular d dimensional sys¬ 
tems, where for any regular simply connected area, the 
boundary area is proportional to , for the CT things 
are more complicated. One expects that a region which 



FIG. I: A Cayley tree with coordination number Z — ?> and 
N — A generations 

is connected to the rest of the graph only at one point (a 
branch of the tree), will show a different EE than a region 
centered at the root of the tree where the boundary area 
grows as ^ (where C is related to the coordination 
number of the tree). In this paper we aim to calculate 
the EE for a clean CT as a stepping stone towards inves¬ 
tigating the EE in a disordered CT. We shall pay special 
attention to the different behavior for different regions 
and the physical meaning of these results. 


II. CAYLEY TREE 

A CT i^ a simple connected undirected graph, with 
no closed loops. A CT of coordination number Z and 
N generations is described as follows. There is a root 
vertex, which we denote (1). The root is linked to Z ver¬ 
texes belonging to the 2nd generation, which we denote 
(2), (3),..., (Z + 1). Each of those is linked to another 
Z — 1 vertexes on the 3rd generation, so that (2) is re¬ 
lated to (Z + 2),..., (Z^ + 1), and so on. Every vertex 
is linked to Z others, except those belonging to the last 
generation, which are linked to only one vertex each. Un¬ 
like the Bethe Lattiee, CT is finite, it’s boundary being 
a non-negligible part of the entire tree: the total number 

of vertexes in the tree is 1 + Z -, while the last 

generation contains Z{Z — 1)^“^ vertexes. Therefore, 
even for Y ^ oc the boundary is import ant, as is well 
known in the context of localization on CT^^^. 
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Another noteworthy property of the CT is it’s lack of 
dimension - since the number of vertexes grows exponen¬ 
tially with N, we cannot define the dimension of this sys¬ 
tem. On the other hand, any two vertexes are connected 
by a single path, what gives the CT a one-dimensional 
character. Both of these features are reflected in the scal¬ 
ing of the EE we obtained. 


III. CT EIGENSTATES 


Eirst we solve the tight-binding Hamiltonian for a clean 
CT. The Hamiltonian can be written, using the notation 
introduced previously, in the form 


z+i N-iz{z-iy ^ z-i 

H = + h.c.)-tY, Y E(A+i-vp.+i+(2-i)(i-i)+(j-i)+ h.c.) 

i=2 g=2 i=l j=l 


r 


( 2 ) 


where is an annihilation (creation) operator on ver¬ 
tex i, and Pg = 2 y is the index of the 

first vertex of the ^’s generation. It will be convenient 
to define a new basis, using the operators de¬ 

fined as follows. We find an orthonormal basis in the 
2nd generation, {^ 2 ’^^^^}^=!, and for each such oper¬ 
ator we define it’s successors in the following gener¬ 


ations as a reflection of in the j’s 

generation. Eor example, if 62 = ^(c 2 — C 3 ), then 




(Eifz+ 2 Ci - Similarly 


^3Z-l 
-'j 


V2^ 


P’"' = 1 

we find a basis in every generation, composed of operators 
other than the successors of the previous generations, and 
define its successors in the next generations. We denote 
each such operator by where j is the generation 

in which this operator acts, T is the generation in which 
it’s first fore-father acts, and u is it’s fore-father’s index 
inside the T’s generation. Among these operators there is 
one important group that originates from the root, which 
we shall denote or 6 ^, {j = 1,..., A^). They are 

symmetric in the sense that they treat each generation 
as one unit, i.e. 




^ =T^ EEp^. ^ Ci for > j > 1 


Cl for j = 1 


yZ(Z-l)i-2 

Using this representation, the Hamiltonian becomes 


( 3 ) 


Thus all terms, except for Hg, are a sum of independent 
one-dimensional lattices. One should note the exponen¬ 
tially increasing number of such lattices. The solution is 
given by the transformation 


4 = 


/ 2 

(j — r + i)k'K 

V - r + 2 ^ 

i=r 

v-r + 2 




(6) 


so that 


H = H,- 

N-i Jv-r+1 

EE E 2 cos 


r =2 k=l 


kn 


N-Ty2 


9r,zyt yUji/ 

4 4 


( 7 ) 


and the spectrum is 


f kn \ 

£fc,r = Icos ■ ( 8 ) 

Eor the symmetric term <§ a similar solution was ob¬ 
tained by Chen et including the boundary condition 
for finite N. 


IV. ENTANGLEMENT ENTROPY 


H = Hs- ^ ^ 

j=2 h 

N-1 j Z{Z-2){Z-iy-^ 


N-1 Z-1 


j=2 v=l 


EE E 

i=3 r=3 


+ h.c.) 


U=1 


( 4 ) 


where the symmetric term 


AT-l 


Hs = -tVzihyp^+ h.c.) - Y 

i=2 

( 5 ) 


In order to calculate the EE we use the relatiorP^ 

5 = ^-A,lnA(-(l-A,)ln(l-A0 (9) 

I 

where {A/} are the eigenvalues of the correlation matrix 
in region A. Eor a one-particle state, the correlation can 
be written in the form C = where is the 

eigenstate in region A. Thus one eigenvector is with 
corresponding eigenvalue A = Since C is Her- 

mitian, any other eigenvector is orthogonal to and 

its corresponding eigenvalue therefore vanishes, so that 
the EE is simply S = —A In A — (1 — A) ln(l — A). We have 
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FIG. 2: A typical eigenstate in a CT with Z = 3, when we 
take one branch (solid line) or a part of it (dashed) to be the 
subsystem. (Color online) 


calculated the EE for several representative one-particle 
states in a CT of Z = 3. A typical state is: 


I^F) = rfT^|0) = 

^ 1 

= V_ - _ 


Sin 


U - l)7r 
N 


E 4 - E 4 

l=Pj+2o-‘^ 


i=Pj 


|0) (10) 


which resides on two of the branches. Eor this state, 
we let one branch be region A (see Eig. |^, and obtain 
A = and S = — In | = ln2, regardless of 

the number of generations. Another possible partition is 
to define region A as a part of one branch (see in Eig. 
[^. The eigenvalue Am for this partition where region A 
is defined as the last M generations of one branch is 


= il + T (i _ + 1)]' 


2N m 


sm 


(f) 


which, in the appropriate limit, is proportional to the 
relative number of generations in region A. 

Eor the more challenging case of a many-particle state, 
the correlation matrix becomes more complicated, so we 
diagonalize it numerically. The many particle state can 
be written as 


1^) = n ‘'40), (11) 

r,u,k 

where T, z/, /c are the indicies of the occupied single par¬ 
ticle states, leading to a correlation matrix 

= {cicm) = (0| n n l<^) • 

r,u,k T'y,k' 

( 12 ) 

Let us first consider the ground state EE of a single 
branch of the CT (see in Eig. for the half-filled case. 
Since a branch is connected to the rest of the system only 


FIG. 3: Ground state EE between one branch and the rest 
for Z = 3 as function of M the number of generations in the 
subsystem. H ere M = iV — 1. The black squares represent an 
even number of generations and the red circles represent odd 
M’s, separated in order to emphasize the alternating term. 
The red lines correspond to 0.15 In (M) + 0.29 for even M’s 
and 0.11 In (M)+0.39 for odd M’s. The blue line in between is 
I In (A^) + 0.24, which corresponds to the asymptotic behavior 
of the EE. 


at one point, one would expect from the area law that 
the EE will be constant. Nevertheless, since the CT has 
no gap, which is similar to the situation in metallic ID 
system, a logarithmic correction is expected. Indeed, we 
find a logarithmic dependence for the CT, but unlike the 
ID system it does not depend on the size of the system, 
but rather on the number of generations N in the CT (or, 
equivalently, on the number of generations M in the sub¬ 
system, where in this case M = N — 1). Specifically, for 
Z = 3, 5^ ~ i \nN+{-y 7 {N), where (-1)^ 7 (N) is 
a decaying alternating term twical to systems with open 
boundary conditions^ (Fig. The pre-factor ^ of the 
logarithmic term is reminiscence of the behavior of an 
open boundary ID system, where if one bisects a s ystem 
of length L in the middle the EE grows as ^ In 
The fact that for the EE on a CT the number of gener¬ 
ations N plays the same role as the length L for a ID 
system stems from the structure of the eigenstates of the 
CT demonstrated in Sec. |III| It is interesting to point 
out that for the CT the ratio of the sites on a branch to 
the overall number of sites is nearly constant (^ 1/^) for 
any generation N larger than 0(1), again similar to the 
situation for the bisected ID system. 

On the other hand, for an excited state we find an 
exponential growth of the EE (Eig. |^, reflecting the 
growth in the number of sites per generation. This is 
expected since for an excited state the EE should obey a 
volume law. Nevertheless, this is in stark difference to the 
one dimensional case for which excited states EE grows 
linearly. Thus, for the EE of a branch of the CT there 
is a huge change between the ground state EE and the 
excited state EE not seen in more standard models. This 
may lead to a huge influence of localization on low-lying 
excitation^^. 

A similar form of dissection is to define a site residing 
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FIG. 4: EE between one branch and the rest of the CT with 
Z = 3 for a highly excited state represented by the black 
squares. The EE grows as 0.39 exp (0.71M) — 0.65 0.39 x 

21 . 02 M _ Q 05 since for an excited state the EE is 

proportional to the volume. 



EIG. 5: Different ways to partition the GT: taking M gener¬ 
ations around the root to be region A (solid circles) or taking 
a branch (big dashed circle) or a sub-branch (small dashed 
circle). 

M generations away from the boundary and all the sites 
connecting it to the boundary as region A, so that A is 
connected to the rest of the tree only at one point (see 
in Fig. just like in the previous case. As expected, 
we get the logarithmic dependence for this case also (Fig. 

i- 

A different obvious way to cut the CT into two different 
regions is to define the root and the sites belonging to the 
first M generations as region A, and the rest of the sites 
belonging to the higher generations as region B (Fig. [^. 
In this case the area of contact between the two regions 
is proportional to Z{Z — 1 )^“^ M > 2, and therefore 
according to the area law we expect the ground state EE 
to be exponentially dependent on M. Indeed, in this case 
the EE behaves as Sa ^ exp (aM) (see Fig. [^. Since 
for such a dissection the volume (number of sites) within 

region A is equal to 1 + ^ {Z — 1)^“^, 

i.e., grows exponentially with M. Also for a highly ex¬ 
cited state we find the EE to follow Sa ^ exp((a'M), 
as expected. Thus, for this form of dissection there is no 
qualitative difference between the EE of the ground state 



EIG. 6 : EE when we cut the GT to sub-branches with M 
generations, where the total number of generations in the tree 
is N = 13. The black squares represent even M’s, with the 
red line corresponding to 0.15 In (M)+0.28, and the red circles 
are for odd M’s, with the line corresponding to 0.086 In (M) + 
0.47. 



M 


EIG. 7: EE between the first M generations around the root 
(including the root itself) and the rest of the GT with Z = 3 
for the ground state (black squares) and for a highly excited 
state (red circles). In both cases the EE grows exponentially 
(the lines, corresponding to 0.55 exp (0.68M) — 0.44 0.55 x 

20.98M _ Q 44 1.04 exp (0.69M) - 1.39 1.04 x 2^ - 1.39, 

respectively), since there is no qualitative difference in the 
growth as function of M between the interface area and the 
volume. 


and the behavior of an excited state as function of M. 


V. DISCUSSION 

In this paper we calculated the EE for the ground state 
and excited states of a clean CT. Utilizing an exact ex¬ 
pression for the single-particle eigenfunctions of the CT 
developed in this paper in order to to obtain a correla¬ 
tion matrix from which the EE of a many-particle state 
can be calculated, it has been demonstrated that the EE 
may show dramatic differences depending on the way the 
CT is partitioned into two regions. If the CT is dissected 
into a branch connected in a single point to the CT, a 
logarithmic dependence on the number of generations in 
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the branch is found, which reflects the one dimensional 
character of the states of the CT. On the other hand, 
excited states on the branch show a completely different 
EE behavior, namely an exponential dependence on the 
number of generations in the branch. This logarithmic 
to exponential change in the entanglement between the 
ground state and the excited state is unique to the en¬ 
tanglement on the CT. On the other hand, the ground 
state EE shows a completely different behavior for a dif¬ 
ferent sort of partition - where the root and the first M 
generations are considered as one region. In this case 
an exponential dependence on M is found, which reflects 
the growth in the number of sites on the surface. In con¬ 
trast to dissecting by branch, here the excited states EE 
(which remains the volume law) does not show a dra¬ 
matic difference when compared to the ground state EE. 

Thus, how one defines the separation into two regions 
has a crucial influence on the behavior of the EE. These 
large fluctuations in the EE could become even stronger 
for other types of complex (random) networks, such as 
small world graphs, Erdos-Renyi graphs and scale free 


network^^. 

These behaviors predict a strong influence of disorder 
on the EE of a CT. It has been showrP^ that the local¬ 
ization length in ID may be inferred from the saturation 
of the EE on the length scale of the localization length. 
Eor the CT one expects that for a region around the root 
and the first M generations where the increase in the 
EE is exponential this tool will be more sensitive than 
for the ID system where the growth is logarithmic. It 
will be also very interesting to compare the effect of the 
localization on the EE of a branch. Will it be similar 
to the ID case? Another question is whether the non 
ergodic phase predicted for the CtP may be detected in 
the behavior of the EE. 
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